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1.  Introduction 


Storm  time  Pc  5  waves  (period  =  150-600  s),  which  perhaps  are  the 
lowest-f requency  and  larges t-amp 1 i tude  among  storm  time  ULF  waves,  were 
recognized  in  the  late  60 's  and  have  since  been  the  subject  of  many 
observational  as  well  as  theoretical  studies.  Early  work  on  this  subject 
can  be  found  in  Brown  et  al .  [1968],  Sonnerup  et  al.[1969],  Hasegawa 
[1969],  and  Lanzerotti  et  al.  [1969]. 

Recent  observations  with  multiple  satellites  and  ground-based 
radar  established  the  variation  of  the  amplitude  and  phase  of  storm  time  Pc 
5  waves  across  the  ambient  magnetic  field.  The  waves  are  observed  to 
propagate  westward  with  respect  to  the  ground  [Walker  et  al .  .  1982;  Allan 
et  al .  ,  1982]  and  geostationary  satellites  [Takahashi  et  al.,  1985a;  Lin 
and  Barfield,  1985],  The  waves  have  a  large  azimuthal  wave  number  (|m|  = 
20-100),  and  their  typical  propagation  velocity  and  azimuthal  wavelength 
are  5-50  km/s  and  1000-10000  km,  respet.  ively.  For  one  case,  the  radial 
extent  of  the  wave  was  observed  to  be  as  large  as  1.7  R^,  [Takahashi  et 
al. ,1985b].  The  oscillation  phase  of  the  wave  field  with  radial  distance  (or 
with  latitude  for  ground-based  observations)  is  nearly  constant  [Walker  et 
al .  ,  1982;  Takahashi  et  al .  ,  1985a]  or  varies  much  slower  than  for  Pc  5 
waves  that  are  externally  driven  field-line  resonance  [Allan  et  al . ,  1982]. 

In  spite  of  its  importance  in  theoretically  discussing  the 
mechanism  for  wave  excitation  [e.g..  Southwood,  1976],  the  field-aligned 


structure  of  storm  time  Pc  5  waves  is  much  less  understood  than  the 
azimuthal  or  radial  structure.  It  is  usually  assumed  that  the  waves  have  a 
standing  (rather  than  propagating)  structure  along  the  field  line,  and 
there  is  some  evidence  to  support  this  assumption  in  the  form  of  the 
quadrature  phase  relation  between  the  mutually  orthogonal  electric  and 
magnetic  field  components  [Walker  et  al .  ,  1983;  Pokhotelov  et  al .  .  1985]. 
The  standing  structure  of  the  wave  has  been  usually  inferred  from 
observations  such  as  the  phase  and  amplitude  relation  among  the  components 
of  the  electric  and  magnetic  fields  of  the  wave  [Walker  et  al.,  1983; 
Pokhotelov  et  al..  1985]  or  particle  flux  modulation  caused  by  the  wave 
[Takahashi  Higbie,  1986].  Only  rare  observations  of  the  nodal 

structure  of  storm  time  Pc  5  waves  have  been  made  with  spacecraft  with 
varying  magnetic  latitude  [Takahashi  et  al .  ,  1985b].  Regarding  the  wave¬ 
length  along  the  field  line,  Nagano  and  Araki  [1983]  argued  that  the  waves 
have  a  significant  amplitude  only  near  the  equator.  Their  argument  was 

based  on  a  statistical  study  of  the  magnetic  field  data  from  GOES  3 

o  o 

(magnetic  latitude  =  5  )  and  GOES  3  (9  ).  Probably  owing  to  insufficient 
spatial  coverage  of  the  data,  these  studies  led  to  various  standing  wave 
models  for  storm  time  Pc  5:  a  symmetric  wave  [Walker  et  al .  ,  1982],  an 

antisymmetric  wave  [Pokhotelov  et  al.,  1985;  Takahashi  and  Higbie,  1986], 
and  a  combination  of  a  symmetric  wave  and  an  antisymmetric  wave  [Walker  et 
al . ,  1983;  Takahashi  et  al.,  1985b]. 

Hereinafter  we  define  the  symmetry  of  a  standing  wave  in  terms  of 
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field  line  displacement.  With  this  definition,  fixed-end,  odd-mode 
(even-mode)  standing  waves,  which  have  an  antinode  (node)  of  field-line 
displacement  at  the  equator,  are  symmetric  (antisymmetric).  Since  the 
transverse  component  of  the  perturbed  magnetic  field  is  proportional  to  the 
tilt  of  the  local  field  line,  a  symmetric  (antisymmetr ic)  wave  has  a  node 
(antinode)  at  the  equator  for  the  transverse  component.  In  contrast,  since 
magnetic  field  compression  is  inversely  proportional  tothe distance  between 
field  lines  (for  one-dimensional  compression),  a  symmetric  (antisymmetric) 
wave  has  am  antinode  (node)  at  the  equator  for  the  compressional  component. 

In  this  report,  we  study  the  field-aligned  structure  of  a  storm 
time  Pc  5  wave  using  magnetic  field  data  from  the  four  satel 1 i tes— SCATHA , 
GOES  2,  GOES  3.  and  GEOS  2  —which  were  located  at  or  near  geostationary 
orbit.  We  have  chosen  a  Pc  5  wave  which  occurred  on  November  14-15,  1979, 
because  the  event  lasted  for  longer  than  50  hours,  and  because  from 
previous  reports  on  the  same  event  [Higbie  et  al . ,  1982;  Takahashi  et  al . , 
1985a],  we  have  some  knowledge  on  the  propagation  characteristics  and  local 
time  extent  of  the  wave. 

In  support  of  a  previous  study  based  on  energetic  particle  data 
[Takahashi  and  Higbie.  1986],  we  obtain  clear  evidence  that  the  wave  had  an 
antisymmetric  standing  structure  along  the  ambient  magnetic  field  line. 
Furthermore,  we  demonstrate  that  the  wave  length  parallel  to  the  field  line 
was  only  a  few  earth  radii,  which  is  much  shorter  than  the  length  of  the 
field  line  itself.  These  observations  should  make  a  considerable  impact  to 


i 
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the  theory  of  storm  time  Pc  5  waves.  The  organization  of  this  report  is  as 
follows.  In  section  2,  experiments  and  data  are  described.  In  Section  3, 
the  data  are  analyzed.  In  Section  4,  a  phenomenological  model  of  a 
standing  wave  is  presented.  In  Section  5,  we  discuss  the  result  of  the  data 
analysis  in  relation  to  the  theories  of  the  propagation  and  generation  of 
hydromagnet ic  waves.  Section  6  is  the  conclusion. 

2.  Experiments  and  Data 

The  magnetic  field  data  have  been  acquired  with  fluxgate 

magnetometers.  The  descriptions  of  the  experiments  have  been  given  by 

Fennell  [1982]  for  SCATHA.  by  Grubb  [1975]  for  GOES  2  and  GOES  3,  and  by 

Knott  [1982]  for  GEOS  2.  We  used  1-min  vector  magnetic  field  data,  which 

were  obtained  by  taking  running  averages  of  the  original 

high-time-resolution  data.  We  also  used  data  from  the  plasma  experiment  on 

board  SCATHA  [Fennell,  1982]  to  determine  ion  pressure. 

We  were  very  fortunate  that  the  four  satellites  were  located  at 

various  magnetic  latitudes.  The  geostationary  satellites  GOES  2,  GOES  3, 

and  GEOS  2  were  located  on  the  geographic  equator  at  the  geographic 
o  o  o 

longitudes  of  107  W,  135  W,  and  15  E,  which  correspond  to  the  magnetic 
0  0  o 

latitudes  of  ~9  N,  ~5  N,  and  ~1  N,  if  we  assume  a  centered  dipole  with  its 

o  o 

north  pole  directed  to  78.3  N  and  291.0  E  in  geographic  coordinates. 
SCATHA  had  an  elliptical  orbit  with  an  orbital  period  of  23.6  hours.  The 
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apogee  and  perigee  were  7.8  R^,  and  5.3  R^.,  respectively,  and  the 

inclination  of  the  orbit  was  7.S  .  During  the  course  of  the  November  14-15 

o  o 

event,  its  magnetic  latitude  varied  between  -7  and  0  . 

3.  Data  Analysis 

Figure  1  summarizes  the  Pc  5  activity  associated  with  the 
geomagnetic  storm  of  November  13-15,  1979.  At  the  top  is  the  hourly  D 
index,  and  at  the  bottom  is  the  local  time  of  the  satellites,  both  plotted 
as  a  function  of  universal  time.  In  the  lower  panel,  the  presence  of  a 
Pc  5  wave  as  observed  by  each  spacecraft  is  indicated  by  the  heavy  portion 
of  the  orbit  traces.  The  wave  lasted  at  least  the  54-hour  period  between 
1800UT  of  November  13  and  2400UT  of  November  15.  This  corresponds  to  the 
recovery  phase  of  the  geomagnetic  storm  shown  in  the  upper  panel.  From 
0700UT  of  November  14  to  0500UT  of  November  15,  the  wave  was  observed  over 
the  entire  dayside;  therefore  it  deserves  the  designation  of  "a  global  Pc  5 
wave"  [Higbie  et  al.,  1982].  We  chose  this  22-hour  interval  for  a  detailed 
study  of  the  spatial  structure  of  the  wave. 


3. 1  Latitudinal  variation  of  phase  and  amplitude 

In  Figure  2  we  show  the  SCATHA  magnetic  field  data  to  demonstrate 
the  continuous  wave  activity.  At  the  top  are  the  1-min  averages  of  the 
magnetic  field  components  By,  B^,  and  B^,  where  e^  (northward)  is 
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Figure  1.  Overall  Illustration  of  Pc  5  Activity  (lower  panel)  and  Its 

Relation  to  the  Ring  Current  Intensity  (upper  panel)  During  the 
Interval  of  13-15  November  1979.  In  the  lower  panel,  the  local 
time  of  the  four  satellites  is  illustrated  as  a  function  of 
universal  time.  The  heavy  portion  of  the  orbit  trace  means  the 
presence  of  a  Pc  5  wave. 
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Figure  2.  SCATHA  Magnetic  Field  Showing  a  Continuous  Wave  Activity  for  th 
12-hr  Period  (upper  panel).  The  lower  three  panels  show  the 
position  of  the  satellite  in  terms  of  magnetic  shell  parameter 
(L),  local  time  (LT),  and  magnetic  latitude  (X). 


antiparal lei  to  the  geomagnetic  dipole,  e^  (eastward)  is  in  the  direction 

of  e„  x  R  ,  where  R  is  the  radial  vector  of  the  satellite  with  respect  to 
H  ~s  ~s 

a  a  /v 

the  center  of  the  earth,  and  e^  (outward)=  e^  x  e^  completes  a  triad.  At 
the  bottom  are  the  L  value,  local  time,  and  the  magnetic  latitude  (A)  of 
the  satellite.  We  used  the  centered  dipole  mentioned  above  to  calculate  L 
and  X. 

The  wave  is  present  continuously  for  the  12-hour  interval  shown. 
It  has  an  average  period  of  600  s,  and  all  the  three  components  show  an 
oscillation  of  the  same  periodicity.  The  H  component,  which  has  the 
largest  amplitude  overall,  demonstrates  the  highly  compressional  nature  of 
the  oscillation.  The  background  field  changes  smoothly  with  the  combined 
effects  of  the  varying  L  and  X. 

For  the  interval  of  1900-2300UT,  we  obtained  =  0.16-0.21  and 
0^  =  0.23-0.30  from  the  SCATHA  ion  data  (17  eV  to  300  keV),  where  0  is  the 
ratio  of  plasma  pressure  to  magnetic  field  pressure,  and  the  subscripts  "II" 
and  "1"  refer  to  the  components  parallel  and  perpendicular  to  the  ambient 
magnetic  field.  These  values  do  not  satisfy  Hasegawa’s  [1969]  condition 
for  drift-mirror  instability 

h 

1  +  0(1  -  ~)  <  0. 
x  P(l 


The  same  interval  is  shown  in  Figure  3  using  an  expanded  time 
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Sequence  of  SCATHA  Magnetic  Field  Data  Presented  in  the  Detrended 
Field-Aligned  Coordinates.  The  position  parameters  given  at  the 
right  of  each  panel  correspond  to  the  center  of  the  2-hr  segment. 
Vertical  dotted  lines  are  drawn  to  indicate  the  maxima  of  bz. 
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scale  and  a  field-aligned  coordinate  system.  In  this  system,  e 

z 

(northward)  is  along  the  ambient  field  averaged  over  a  1-hour  interval 

A  ^  A 

about  each  1-min  data  point,  e  (eastward)  is  parallel  to  e z  x  £s>  and 

A  A 

(outward)  =  e  x  e  completes  the  triad.  For  the  z  component,  the  trend 

y  ^ 

has  been  removed  by  running  average.  In  the  text,  the  perturbation  in  the 

field  components  shall  be  denoted  as  b  .  b  ,  or  b  hereinafter.  We  have 

x  y  z 

organized  the  five  2-hour  segments  according  to  the  geomagnetic  latitude  of 

the  satellite,  so  the  time  runs  from  bottom  to  top.  The  L  value,  local 

time,  and  X  are  given  for  the  center  of  each  2-hour  segment.  In  Figure  4 

we  have  added  the  magnetic  field  data  from  the  three  geostationary 

satellites  located  at  fixed  magnetic  latitudes.  The  format  is  the  same  as 

in  Figure  3.  All  the  data  were  taken  around  noon. 

The  key  features  to  be  noticed  in  Figure  3  and  Figure  4  are  the 

phase  and  amplitude  relations  among  the  field  components.  For  example,  by 

leads  b  by  90  at  X  =  -6.5°  and  -5.7°.  b  has  a  small  amplitude  at  X  = 
x  y 

o  o 

-4.1  so  that  the  relative  phase  cannot  be  defined,  and  at  X  =  -2.6  and 

O  0  0 

1.2  .  b^  leads  b^  by  90  .  At  X  =  4.7  ,b  again  has  a  very  small  amplitude, 
o 

At  X  =  9.2  ,  b^  has  a  small  amplitude,  making  it  difficult  to  establish  the 
phase  lag.  Amplitude  and  phase  variations  are  also  observed  for  b^. 

Although  the  wave  event  lasted  very  long,  it  is  still  possible 
that  the  amplitude  changed  with  universal  time,  local  time,  or  L  value,  as 
well  as  with  magnetic  latitude.  Thus,  use  of  raw  amplitude  makes  the 
modeling  of  the  field-aligned  structure  somewhat  difficult.  For  this 
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reason  we  introduced  a  normalized  amplitude 
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to  describe  the  variation  of  relative  amplitudes  among  the  field  components 


with  magnetic  latitude.  Here  a  =  x,  y,  or  z,  and  <b^>  is  the  mean-square 


amplitude  of  the  magnetic  field  component  obtained  by  integrating  its  power 
spectral  density  over  the  frequency  range  occupied  by  the  wave.  The 
spectral  analysis  was  done  on  each  of  the  2-hour  segments  shown  in  Figures 
3  and  4.  In  the  same  spectral  analysis,  the  relative  phases  among  the 


magnetic  field  components  were  also  obtained.  We  denote  by  <p^  the 


relative  phase  of  component  a  with  respect  to  component  0. 

The  normalized  amplitude  and  the  relative  phase  are  illustrated 
in  Figure  5  and  Figure  6.  respectively,  as  a  function  of  magnetic  latitude. 
The  circles  are  the  observational  data  points,  and  the  curves  are  given  by  a 
model  standing  wave  described  below.  In  Figure  5,  it  is  clear  that  the 
transverse  components  have  a  maximum  amplitude  near  the  magnetic  equator 
and  a  minimum  several  degrees  off  the  equator.  In  contrast,  the 
compressional  component  has  a  minimum  amplitude  near  the  equator.  In 


Figure  6,  it  is  seen  that  the  relative  phase  $  or  *  is  either  -90  or 
®  r  xy  yz 


O  0  0 

90  ,  and  a  switching  occurs  somewhere  between  X  =  -5  and  -3  .  The  relative 


O  0 

phase  $  is  either  0  or  180  ,  and  a  switching  occurs  somewhere  between  X  = 


U 


a 


i: 


a 


'•j 
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Normalized  amplitud 


We  can  argue  that  only  has  a  clear  indication  of  a  node  within 
the  latitude  range  covered  by  the  SCATHA  data.  Figure  7  illustrates  the 
latitude  dependence  of  the  normalized  amplitude  of  the  transverse  magnetic 
field  components  and  the  phase  lag  between  them.  Note  that  we  used  only 
the  transverse  components  for  the  normalization  in  order  to  emphasize  the 
difference  between  the  x  and  y  components.  That  is.  we  have  defined 


<b*>  1 

( - 2 - )  1/2 

v  p  O  ' 

<b  p  +  <b> 

x  y 


a  =  x  or  y. 


2  2 

Also,  we  used  the  mean  variance  <bx>  or  <b^>  calculated  on  a  cycle-by-cycle 

basis  as  an  indicator  of  the  amplitude,  where  the  wave  cycle  was  defined  as 

the  interval  between  two  successive  peaks  of  b^.  This  method  allows  us  to 

determine  the  latitude  structure  of  the  wave  in  a  fine  scale  and  also  to 

evaluate  the  degree  of  data  scattering  at  a  given  magnetic  latitude.  The 

o 

middle  panel  indicates  that  by  has  a  relative  minimum  at  X  ~  4  ,  and  the 

lower  panel  indicates  that  a  change  in  relative  phase  occurred  somewhere 
o  o 

between  X  =  -4.2  and  X  =  -3.3  .  These  observations  can  be  taken  as 

o 

resulting  from  a  node  of  b  located  near  -4  .  The  value  of  $  is  constant 

y  xy 

o  o  o  _  o 

between  X  =  -6.6  and  X  =  -4.2  ,  and  also  between  X  =  -3.3  and  X  =  -1.2  , 

from  which  we  conclude  that  there  is  no  node  for  b  or  b  in  these  latitude 

x  y 

ranges.  This  conclusion  is  further  confirmed  when  we  compare  the  amplitude 
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Figure  7.  Cycle-by-Cycle  Amplitude  and  Phase  of  the  Transverse  Magnetic 
Field  Components.  The  vertical  lines  indicate  a  possible 
latitudinal  range  of  the  nodes  of  a  standing  wave.  The  relative 
phase  4>  was  determined  by  eyeball  examination  of  the  time  plots 
of  bx  ar$  b  .  Whenever  b  (b  )  led  b  (bx)  by  a  phase  difference 
smaller  than  180°,  the  value  of  was  set  to  90°  (-90°). 


plots  between  the  two  components.  The  minimum  value  of  b  ,  corresponding 
to  a  b  node,  is  ~0.2.  not  exactly  0.  This  can  be  regarded  as  the  result 

y 

of  noise  added  to  the  wave.  Assuming  that  the  noise  has  the  same  amplitude 

for  b  ,  we  would  expect  that  b  would  also  reach  ~0.2  at  its  node.  All  the 
x  x 

data  points  for  the  component  shown  in  the  upper  panel,  however,  are 
greater  than  0.3;  therefore  we  conclude  that  no  node  of  b^  exists  between  A 
=  -6.6°  and  A  =  -1.2°. 

From  an  inspection  of  the  GOES  2  data,  on  the  other  hand,  we  can 

O 

conclude  that  a  node  of  b  is  located  near  the  GOES  2  latitude  of  9.2  .  In 

x 

Figure  2,  the  overall  amplitude  of  the  x  component  at  GOES  2  is  almost 

null.  As  an  exception,  the  three  wave  cycles  between  1825  and  1845  UT 

indicate  a  small  but  finite  amplitude  of  b^,  and  during  this  period  bx 
o 

leads  by  by  ~90,  which  is  the  same  as  the  observation  near  the  equator  made 

by  SCATHA  or  GEOS  2.  This  phase  relation  can  be  explained  by  assuming  that 

a  node  was  located  between  the  equator  and  GOES  2  for  both  b  and  b  .  since 

x  y 

O 

♦  changes  by  180  each  time  a  node  is  crossed.  However,  we  maintain  the 
xy 

view  that  b  and  b  have  different  node  latitudes,  as  was  the  case  with  the 
x  y 

SCATHA  observation.  An  inspection  of  the  ratio  between  b  and  b  in  the 

x  y 

upper  two  panels  of  Figure  3  leads  to  the  conclusion  that  b^  has  a  node 
near  GOES  3  (5  ) .  while  b^  has  a  node  near  GOES  2  (9°). 


3.2  Harmonics 


An  interesting  wave  form  is  observed  in  b  and  b  .  An  example  is 

y  ^ 


k 
k 

I 

i 

1 

k 

I 

1  seen  in  the  bz  trace  of  the  lower  panel  (GEOS  2)  of  Figure  4.  In  this  case 

!  b  has  an  amplitude  smaller  than  the  transverse  components,  but  we  can  see 

I 

|  that  it  oscillates  at  twice  the  frequency  of  b  and  b  .  The  b  wave  form 

,  x  y  z 

1  is  not  a  simple  sinusoid:  two  consecutive  peaks  do  not  have  the  same 

!  height  except  near  1050UT. 

Coleman  [1970],  who  analyzed  ATS  1  magnetic  field  data,  was  the 
first  to  report  on  this  type  of  wave  form.  He  called  it  the  "rectified” 
version  of  b^.  Kokubun  [1985]  found  the  same  phenomenon  in  GOES  2  and  GOES 
3  data  and  called  it  a  "harmonic."  Kokubun  speculated  that  this  phenomenon 
is  localized  in  magnetic  latitude,  because  the  harmonic  was  more  often 
observed  at  GOES  3  (X  ~  5°)  than  at  GOES  2  (X  ~  9°). 

Our  data  set  from  the  four  satellites  clearly  demonstrates  that 
Kokubun 's  speculation  is  correct:  the  harmonic  bz  oscillation  occurred 

O 

only  at  X  ~  1  in  our  case.  The  persistence  of  this  phenomenon  and  its 
stability  in  latitudinal  localization  was  confirmed  by  scanning  the  data 
for  0600-2400UT  of  November  15,  the  interval  following  the  one  analyzed 

O 

above.  Only  GEOS  2  (X  ~  1  ),  among  the  four  satellites,  observed  the 
harmonic . 

A  careful  inspection  of  the  GOES  3  data  in  Figure  4  (its  expanded 


version  is  shown  in  Figure  8)  reveals  that  a  similar  wave  form  occurs  at  X 


~  5  in  by.  Although  the  wave  form  is  essentially  the  same  as  for  bz>  in 
this  case  it  is  not  appropriate  to  call  b^  the  rectified  version  of  any 
other  components  because  the  peaks  occur  at  epochs  when  b^  and  bz  are  close 


4.  Phenomenological  Model 


A  simple  phenomenological  standing  wave  model  can  be  constructed 
to  describe  the  observed  latitude  dependence  of  the  amplitude  and  phase  of 
the  wave.  The  model  we  adopted  and  show  in  Figures  5  and  6  has  the 
following  form: 


bx=  <3xC0S{kz^'ZO*}e^kyy*Ut^ 

by=  iPycos{kz(z-zo)}el{kyy 
b  =  -C  sin{k'(z-z  )}el{kyy_a;t) 


where  (a  =  x.  y,  or  z),  kz>  kz>  and  u  are  real  and  positive.  The  phase 
factor  e^ky^  <l)t^  takes  into  account  the  observed  azimuthal  propagation  of 
the  wave.  According  to  Takahashi  et  al .  [1985a],  who  studied  the  same  wave 
"27r  -1 

event,  k^~  q  g  (R^.  ),  where  the  minus  sign  means  westward  propagation. 

The  relative  phase  among  the  components  at  any  latitude  (i.e.,  z 
coordinate)  can  be  derived  by  taking  the  real  part  of  (la)-(lc). 

The  parameters  in  (la)-(lc)  are  determined  by  trial  and  error.  A 

2 '.IT 

reasonable  fit  to  the  data  was  obtained  with  0:0:6=  4:3:8,  k  =  ~ z 

x  y  z  z  3.0 

- 1  ^  -2 

(Rj.  ).  kz  =  ^-3  (^£  )•  and  z0=  0  085  R^,  where  zq  accounts  for  the 


apparent  offset  of  the  wave  node  from  the  nominal  magnetic  equator.  The 

fitting  of  the  model  was  done  using  the  2-hour  data  points.  Although  the 

model  has  nodes  for  b  and  b  at  latitudes  slightly  different  than  those 

x  y 

discussed  in  section  3  using  cycle-by-cycle  wave  analysis,  the  important 
fact  that  the  different  components  can  have  different  node  latitudes  is 
well  taken  into  account.  In  practice,  zq  is  much  smaller  than  the  parallel 
wavelength,  so  we  can  regard  that  the  node  (anitnode)  of  bz  (b^  and  b^) 
lies  at  the  equator.  This  means  that  the  wave  has  an  antisymmetric 
standing  wave  structure. 

Because  we  are  trying  to  model  the  wave  over  the  latitudinal 

O  O 

range  from  -7  to  9,  we  can  use  an  approximation  z  ~  LR^X  to  relate  the 
magnetic  latitude  to  distance  along  the  field  line.  We  assumed  in  Figure  5 
and  Figure  6  that  the  relative  phase  and  the  normalized  amplitude  are  a 
function  of  magnetic  latitude  only.  Thus  the  corresponding  wavelength 


depends  on  the  L  value.  The  values  for  k  ,  k',  and  z  used  above 

z  z  o 


correspond  to  the  geostationary  distance  L  ~  6.6.  We  should  also  emphasize 
that  the  model  (la)-(lc)  probably  applies  only  to  the  latitudinal  range  of 
X  ~  -10°  to  X  ~  10°,  in  which  geostationary  satellites  stay.  At  present  we 
have  no  knowledge  how  the  wave  structure  might  change  for  |x|  £  10°. 
However,  it  is  obvious  because  of  the  inhomogeneity  of  the  ambient  magnetic 
field,  we  cannot  simply  assume  that  the  wave  field  changes  with  the  same 
wavelength  all  the  way  to  the  ionosphere. 
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5.  Discussion 


The  observations  presented  above  and  summarized  in  terms  of  a 
simple  standing  wave  model  have  the  following  physical  implications. 


5.1.  Short  parallel  wavel eng th 

We  have  clearly  observed  nodes  for  the  transverse  components  at 
several  degrees  off  the  dipole  equator.  This  means  that  the  parallel  wave 
length  is  only  a  few  earth  radii  for  the  components.  We  have  observed  only 
one  node  for  the  compressional  component , but  this  is  because  of  the  limited 
latitude  coverage  of  the  satellites.  It  is  very  likely  that  the 
compressional  component  also  has  a  parallel  wavelength  of  a  few  earth 
radii.  Such  a  short  parallel  wavelength  is  not  consistent  with  the  model 
in  which  a  compressional  wave  (drift-mirror  mode)  is  coupled  to  the 
fundamental  mode  of  the  guided  poloidal  wave  (shear  Alfven  mode)  as  was 
suggested  by  Walker  et  al .  [1982].  A  fundamental  mode  has  half  a  wave¬ 

length  along  the  geomagnetic  field  line.whicn  means  a  wavelength  of  30  R^. 
or  so  at  L  ~  6.6.  In  fact,  with  the  wave  frequency  of  0.002mHz  and  the 
parallel  wavelength  (~3  R^)  obtained  above,  any  shear  Alfven  resonance  u 

t 

-  kzVA  =  0  is  not  likely,  where  is  the  Alfven  velocity.  From  the 

frequencies  of  multiharmonic  standing  Alfven  waves,  the  equatorial  plasma 

mass  density  during  the  November  14  event  has  been  estimated  to  be  25 
3 

amu/cm  [Takahashi  et  al.,  1985a].  For  this  density,  the  Alfven  wave  must 
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The  two  transverse  components,  b  and  b  ,  have  nodes  at  different 

x  y 

magnetic  latitudes.  The  same  phenomenon  was  observed  by  the  ISEE-1  and  -2 
spacecraft  by  Takahashi  et  al.  [1985b],  They  explained  the  observation  by 
invoking  a  symmetric  wave  for  b  and  an  symmetric  wave  for  b  .  With  our 

y  5c 

present  observation,  we  conclude  that  it  is  not  necessary  to  assume  the 
symmetric  and  antisymmetric  combination.  All  the  components  oscillate  in 
an  antisymmetric  mode,  but  the  wavelength  is  component  dependent. 

Southwood  and  Kivelson  [1986]  theoretically  demonstrated  that  a 
compressional  wave  and  a  transverse  wave  can  have  different  amplitude 
variations  along  the  field  line,  provided  that  the  plasma  density  has  a 
gradient  along  the  field  line.  Although  in  their  argument  a  box-model  of 
the  magnetosphere  and  a  cold  plasma  were  assumed,  the  essential  ingredient 
of  the  theory,  the  field  aligned  density  gradient,  should  be  generally 
present  in  the  real  magnetosphere:  therefore,  their  theory  may  apply  to  our 
result . 


5.3.  Polarization 

The  magnetic  field  perturbation  of  the  wave  has  components  in  the 


x,  y,  and  z  directions.  The  relative  phase  among  the  magnetic  field 

o  o 

components,  as  illustrated  in  Figure  6,  is  either  0,  ±90  ,  or  180  and 

O 

switches  by  180  as  the  observing  satellites  move  across  the  nodes  of  the 
wave.  Here  we  demonstrate  that  these  polarization  characteristics  can  be 
explained  in  terms  of  wave  propagation  in  an  inhomogeneous  medium.  What 
follows  is  a  slight  modification  to  the  discussion  given  by  Walker  et  al . 


[1982], 


In  general,  the  perturbation  electric  field  E  for  a  wave  having  a 


wave  vector  k  satisfies 


k  (k*E)  -  k2!  +  (I  +  e)-E  =  0 


where  e  is  the  dielectric  tensor.  Walker  et  al .  [1982]  derived  the 

expression  of  e  for  a  wave  propagating  perpendicular  to  the  gradient  of  the 

magnetic  field  strength  in  a  plasma  which  is  characterized  by  ion  guiding 

center  drift.  Using  their  result,  we  find  that  the  transverse  electric 

field  components  E  and  E  are  related  as 
x  y 


■[(P„  -  p±  -  2)k*  ♦  ^  JEy  ♦  ikyK(P||  *  p±)Ex  =  0. 

A 


We  have  neglected  electron  pressure,  which  is  small  in  the  storm  time  ring 
current  region.  The  parameter  K,  which  is  real  and  positive,  is  equivalent 


* 


to  the  inverse  gradient  scale  length  of  the  magnetic  field  intensity. 

Since  a  local  approximation  has  been  made,  equation  (2)  is  valid  for  the 

case  |ky|  >>  K,  6^,  where  is  the  radial  scale  length  of  the  wave 

field.  For  a  wave  field  which  varies  as  exp[i(k^y  +  k^z  -ut)],  we  obtain 

from  Faraday's  law 
a) 

-k  E  =  ~  b  (3a) 

z  y  c  x  ' 


k  E  =  “  b 
z  x  c  y 


-k  E  =  -  b 
y  x  c  z 


where  3Ey/3x  =  0  is  used  in  the  spirit  of  local  approximation.  Also  we 

assumed  E  =0. 
z 

In  the  manner  described  by  Walker  et  al.  [1982]  we  can  construct 
a  standing  wave  by  the  superposition  of  two  propagating  waves.  Unlike 
Walker  et  al.  [1982],  who  assumed  a  symmetric  mode,  we  shall  assume  an 
antisymmetric  mode  for  the  standing  wave.  In  this  model,  we  take  that  the 
radial  magnetic  field  perturbation  consists  of  the  following  two  wave 
components 


where  A  is  a  constant  representing  the  amplitude  of  the  waves.  Omitting 
the  common  phase  factor  e  Ky- 

,  the  perturbations  in  other  field 
components  corresponding  to  the  above  waves  can  be  obtained  from  (2)  and 
(3a)-(3b)  as 


b*. 
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•  ,2 

iak 

Z  ± 

inrbx 
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iak 

b*  =  i 
z  K  x 


where  we  have  def ined 


2u 

-«P,I  -  h  -  2>  *  ^2  1 

A  z 


a  = 


(*n  +  *i> 


(4a) 


(4b) 


Then,  the  standing  wave  structure  for  each  component  becomes 


b  =  b+  +  b  =  Acos(k  z) 
X  X  X  z 


(5a) 
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+  b 
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Acos(kzz) 
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=  b+  +  b" 
z  z 


ak 

z 

~  Asin(kzz)  . 


(5b) 


(5c) 
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From  the  discussion  presented  above,  k  >>  ~T  holds  for  our  observation, 

z  vr*£ 

A 

and  from  SCATHA  plasma  experiment,  we  have  p^  =  0.23-0.30  and  P^  = 
0.16-0.21.  Those  taken  into  account,  we  find  a  >  0.  Then,  taking  k^  <  0 
to  be  consistent  with  the  observed  westward  propagation,  we  find  that 
(5a)-(5c)  have  the  same  form  as  (la)-(lc)  except  for  the  different  parallel 
wavelengths  for  different  components  and  a  small  shift  of  the  node  from 
the  equator.  Thus  if  we  ignore  these  minor  points,  the  observed  phase 
relations  among  the  components  can  be  explained  by  the  technique  of  Walker 
et  al.  [1982], 


Walker  et  al .  associated  b  with  a  shear  Alfven  wave  and  b  and 


x  y 

bz  with  a  drift-mirror  wave.  They  argued  that  the  drift-mirror  wave 
propagates  at  the  guiding  center  drift  velocity  of  ions  instead  of  the 
diamagnetic  drift  velocity  as  was  originally  discussed  by  Hasegawa  [1969]. 
This  argument  has  some  observational  support:  our  wave  propagated  westward 
at  a  velocity  of  ~10  km/s,  the  guiding  center  drift  velocity  expected  for 
ring  current  ions  with  energy  of  ~10  keV.  However,  the  drift-mirror 
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instability  condition  was  not  satisfied,  at  least  at  the  SCATHA  location. 

Walker  et  al .  further  argued  that  the  wave  is  excited  when  the  coupling 

condition  u  =  k  V,  =  k  V.  between  the  drift-mirror  wave  and  the  Alfven  wave 
yd  z  A 

is  satisfied,  where  is  the  guiding  center  drift  velocity  of  ions.  As 
discussed  above,  this  prediction  does  not  match  our  observation  either,  as 
discussed  above.  It  seems  that  Walker  et  al.'s  theory  has  a  limited 
success  here. 

5.4.  Antisymmetric  standing  wave  structure 

Our  Pc  5  wave  clearly  had  an  antisymmetric  standing  wave 

structure.  At  present  we  do  not  know  why  this  is  the  case.  Drift  mirror 

instability  of  Hasegawa  [1969]  or  Walker  et  al .  [1982]  does  not  have  a 

mechanism  to  excite  an  antisymmetric  wave  preferentially,  because  the  slab 

geometry  which  they  assumed  does  not  take  into  account  the  presence  of 

plasma  and  field  inhomogeneity  in  the  direction  of  the  ambient  magnetic 

field  or  the  presence  of  bouncing  particles. 

In  view  of  bounce-drift  resonance  excitation  of  ULF  waves  in  the 

earth’s  dipole  field  [Southwood,  1976],  there  is  an  intimate  relation 

between  the  field-aligned  structure  of  the  wave  and  the  type  of  resonance. 

On  one  hand  there  is  drift  resonance  u  -  mo,  =  0,  which  could  lead  to 

d 

generation  of  symmetric  standing  waves.  Here,  m  is  the  azimuthal  wave 
number,  and  is  the  drift  frequency  of  a  given  ion  species.  This 
preference  for  symmetric  waves  does  not  allow  us  to  apply  the  resonance  to 
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our  observation,  although  the  wave  angular  velocity  u/m  =  predicted  by 
the  resonance  condition  could  explain  why  the  wave  propagates  westward  at  a 
velocity  ( — 10  km/s)  comparable  to  the  guiding  center  drift  velocity  of  the 
ring  current  ions  (~10  keV). 

On  the  other  hand,  the  drift-bounce  resonance  o>  -  mo,  =  ±o,  is 

d  b 

related  to  generation  of  ant isymme trie  waves,  where  is  the  bounce 

frequency.  Because  of  the  added  term  on  the  right-hand  side  of  the 

equation,  this  resonance  condition  does  not  always  imply  westward 

propagation  (i.e.,  o/m  <  0  ).  However,  we  could  estimate  the  energy  of 

protons  which  is  consistent  with  the  observed  values  of  u  ~  2  mHz  and 

m  ~  -60  [Takahashi  et  al..  1985a].  Using  the  formula  of  Hamlin  et  al. 

[1961]  for  particles  mirroring  near  the  equator,  we  find  that  protons  of 

~0.3  keV  or  ~400  keV  can  satisfy  the  resonance  condition  u  -  mu,  =  ±u.\  . 

a  b 

Such  energy  is  outside  the  typical  energy  range  of  ring  current  protons; 


hence  this  resonance  mechanism  condradicts  the  view  that  the  ultimate 
energy  source  for  storm  time  Pc  5  waves  is  the  ring  current  ions. 


5.5  Harmonics 

If  a  bounded  plasma  has  an  MHD  eigenmode  (u,  k) ,  its  harmonic 
(2u,  2k)  is  also  an  eigenmode.  So,  a  simple  way  to  explain  the  observed 
harmonic  is  to  assume  a  smal  1-ampl i tude  wave  with  2u  and  2k.  For  example, 
we  may  have 
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superposed  on  the  wave  described  by  (la)-(lc),  where  6  <<  P  and  6  <<  P  . 

y  y  z  z 


We  assumed  that  each  component  of  (6a)  and  (6b)  has  the  symmetry  opposi te  to 
that  of  its  counterpart  in  (lb)  and  (lc).  This  is  necessary  in  order  to 
have  a  finite-amplitude  harmonic  at  the  nodes  of  the  "fundamental"  wave 
described  by  (la)-(lc),  but  we  cannot  offer  an  explanation  why  this  is  so. 

Such  a  description  is  highly  phenomenological,  and  there  could  be 
other  ways  than  (6a)  and  (6b)  to  describe  the  observed  harmonics  [Higuchi 
et  al..  1986;  Takahashi  et  al .  ,  1987].  The  questions  of  whether  the 
harmonic  is  a  general  property  of  storm  time  Pc  5  waves  and  whether  it  has 
any  fundamental  importance  to  the  wave  excitation  mechanism  need  to  be 
answered  by  further  observational  and  theoretical  work. 


6.  Conclusion 


In  conclusion  we  have  presented  convincing  evidence  that  the  Pc  5 
wave  of  November  M-15,  1979,  had  an  antisymmetric  standing  wave  structure 
with  a  parallel  wavelength  of  a  few  earth  radii.  Although  the  event  was 
unusual  in  that  it  lasted  for  such  a  long  period  as  ~50  hours,  we  speculate 
that  the  structure  could  be  common  to  other  storm  time  Pc  5  waves. 
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Previously  proposed  theories  on  storm  time  Pc  5  waves  did  not 
take  into  account  an  antisymmetric  standing  wave  structure  or  a  short  wave 
length  along  the  ambient  field  lines.  Naturally  they  fail  to  give  an 
explanation  for  our  present  observation.  A  new  theoretical  development 
seems  to  be  required  for  the  excitation  mechanism  of  storm  time  Pc  5  waves. 
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LABORATORY  OPERATIONS 


The  Aerospace  Corporation  functions  as  an  “architect-engineer”  for 
national  security  projects,  specializing  in  advanced  military  space  systems. 
Providing  research  support,  the  corporation's  Laboratory  Operations  conducts 
experimental  and  theoretical  investigations  that  focus  on  the  application  of 
scientific  and  technical  advances  to  such  systems.  Vital  to  the  success  of 
these  investigations  is  the  technical  staff's  wide-ranging  expertise  and  its 
ability  to  stay  current  with  new  developments.  This  expertise  is  enhanced  by 
a  research  program  aimed  at  dealing  with  the  many  problems  associated  with 
rapidly  evolving  space  systems.  Contributing  their  capabilities  to  the 
research  effort  are  these  individual  laboratories: 


Aerophysics  Laboratory:  Launch  vehicle  and  reentry  fluid  mechanics,  heat 
transfer  and  flight  dynamics;  chemical  and  electric  propulsion,  propellant 
chemistry,  chemical  dynamics,  environmental  chemistry,  trace  detection; 
spacecraft  structural  mechanics,  contamination,  thermal  and  structural 
control;  high  temperature  thermoraechanics ,  gas  kinetics  and  radiation;  cw  and 
pulsed  chemical  and  exciraer  laser  development  including  chemical  kinetics, 
spectroscopy,  optical  resonators,  beam  control,  atmospheric  propagation,  laser 
effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  chemical  reactions, 
atmospheric  optics,  light  scattering,  state-specific  chemical  reactions  and 
radiative  signatures  of  missile  plumes,  sensor  out -of -f ield-of -view  rejection, 
applied  laser  spectroscopy ,  laser  chemistry,  ^aser  optoelectronics,  solar  cell 
physics,  battery  electrochemistry,  space  vacuum  and  radiation  effects  on 
materials,  lubrication  and  surface  phenomena,  thermionic  emission,  photo¬ 
sensitive  materials  and  detectors,  atomic  frequency  standards,  and 
environmental  chemistry. 

Computer  Science  Laboratory:  Program  verification,  program  translation, 
performance-sensitive  system  design,  distributed  architectures  for  spaceborne 
computers,  fault-tolerant  computer  systems,  artificial  intelligence,  micro¬ 
electronics  applications,  communication  protocols,  and  computer  security. 


Electronics  Research  Laboratory:  Microelectronics,  solid-^tate  device 
physics,  compound  semiconductors,  radiation  hardening;  electro-optics,  quantum 
electronics,  solid-state  lasers,  optical  propagation  and  communicat tons ; 
microwave  semiconductor  devices,  microwave/oilliraeter  wave  measurements, 
diagnostics  and  radiometry,  microwave/millimeter  wave  thermionic  devices; 
atomic  time  and  frequency  standards;  antennas,  rf  systems,  electromagnetic 
propagation  phenomena,  space  communication  systems. 


Materials  Sciences  Laboratory:  Development  of  new  materials:  metals, 
alloys,  ceramics,  polymers  and  their  composites,  and  new  forms  of  carbon;  non¬ 
destructive  evaluation,  component  failure  analysis  and  reliability;  fracture 
mechanics  and  stress  corrosion ;  analysis  and  evaluation  of  materials  at 
cryogenic  and  elevated  temperatures  as  well  as  in  space  and  enemy-induced 
environments . 


Space  Sciences  Laboratory:  Magnetospheric,  auroral  and  cosmic  ray 
physics,  wave-particle  interactions,  magnetospheric  plasma  waves;  atmospheric 
and  ionospheric  physics,  density  and  composition  of  the  upper  atmosphere, 
remote  sensing  using  atmospheric  radiation;  solar  physics,  infrared  astronomy. 
Infrared  signature  analysis;  effects  of  solar  activity,  magnetic  storms  and 
nuclear  explosions  on  the  earth's  atmosphere,  ionosphere  and  magnetosphere; 
effects  of  electromagnetic  and  particulate  radiations  on  space  systems;  space 
instrumentation. 
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